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Abstract. Thin films of high temperature garnet materials such as yttrium aluminum 
garnet (YAG) doped with rare earths are currently being investigated as selective emitters. 
This paper presents a radiative transfer analysis of the thin film emitter. From this 
analysis the emitter efficiency and power density are calculated. Results based on 
measured extinction coefficients for erbium-YAG and bolmium-YAG are presented. 

These results indicate that emitter efficiencies of 50% and power densities of several 
watts/cm2 are attainable at moderate temperatures (< 1750K). 


INTRODUCTION 

There are two emitter designs that can be used in a thermophotovoltaic (TPV) 
energy conversion system. In a selective emitter system most of the emitted 
photons have energy greater than the bandgap energy of the photovoltaic (PV) cells. 
In a thermal emitter system a narrow bandpass filter is used to shape the emission 
spectrum. The ideal filter allows all photons within the band to pass from emitter to 
the PV cells and reflects all photons outside the band back to the emitter. Both the 
selective emitter system and the filter system are feasible. In this paper we present 
the radiative transfer theory necessary to calculate the spectral emittance of a thin 
film selective emitter. From the spectral emittance, the emitter efficiency is 
calculated. 

The solid state selective emitters of most interest are compounds containing the 
rare earths. For doubly and triply charged ions of these elements in crystals, the 
orbits of the valence 4f electrons, which account for most of the emission and 
absorption, lie inside the 5s and 5p electron orbits. As a result, the rare earth ions 
in the solidstate have radiative characteristics much like they would have if they 
were isolated. They emit in narrow bands rather than in a continuum as do most 
solids. The 5s and 5p electrons "shield" the 4f valence electrons from the 
surrounding ions in the crystal. The spectra of these rare earth ions in crystals have 
been extensively studied. Most of this work is summarized in the text of Dieke(l ). 

Early spectral emittance work(2) on rare earth oxides showed strong emission 
bands. However, the emittance for photon energies below the bandgap for PV 
materials was also significant. As a result, the efficiency of these emitters was low. 
In the last few years, however, large improvement in rare earth emitters has 
occurred. The Welsbach mantle type emitters of Nelson and Parent(3,4) have a 
single strong emission band with greatly reduced emission outside the emission 
band. The very small characteristic dimension (5-10 pm) of the emitter results in 
low emittance outside the emission band. A new thin film rare earth-YAG (ytterium 
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aluminum garnet) emitter(5) shows great promise. The emittance theory presented 
in this paper applies to any thin film emitter. However, the emitter efficiency results 
are calculated specifically for the rare earths since a three band model is used in 
those calculations. In this model the spectrum is split into three wavelength 
regions; below emission band radiation which is described by a constant extinction 

coefficient, a u , within the emission band where a x depends on wavelength, X, and 
above the emission band which is described by a constant extinction coefficient, cc t 
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Substrate emittance = 

y Ztzzzzzzzzzzzzzzzzzzzzzzzzzzi 


= reflectivity of film-substrate interface 
Px* = reflectivity of film-vacuum interface 
K^d = a^d - optical depth of film 


©ji> = black body emissive power 
£^ = spectra! emittance of substrate 
a* = extinction coefficient 


FIGURE 1. Schematic of Emitting Rim 


RADIATIVE TRANSFER THEORY 

The thin film is approximated by the one dimensional model shown in figure 1 
where thermal energy enters through the film substrate. All the film variables such 
as radiation intensity, i x , are assumed to vary only in the x direction. This is an 

excellent approximation for thin films where the film dimensions perpendicular to 
the x direction are much larger than the film thickness, d. The steady state one 

dimensional energy equation 6 for the solid state film is 
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( 1 ) 


-ilk — 

dx[ 1 dx 



where T e is the film temperature, k t is the film thermal conductivity and Q is the 
total radiation energy flux. 



( 2 ) 

(3) 


The spectral intensity, i x , in terms of the wavelength, X, is determined by the 
radiative transfer equation(6). 


1 di ^ • _c 
a x dr lx x 


(4) 


which gives the variation of the intensity in some direction, r, in terms of the so- 
called source function^ and the extinction coefficient, a x . 


«x =a x +cf x. (5) 

where a x is the absorption coefficient and a x is the scattering coefficient. The 
source function includes both emission and scattering and depends on the scattering 
albedo defined as 


^ a x + a x a x (6) 

In this analysis we assume that T E is constant throughout the film although 

experimental results for the rare earth-YAG thin film emitter(5) show a temperature 
gradient across the emitter (200 - 300°K). An analysis that includes temperature 
variation is reported in another paper 8 at this meeting. With T E = constant, the 
limiting solution to eq. (1) is q = constant and all the thermal input leaves the film as 
radiation. The details of the solution for the source function. 
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S^, and intensity, i^, are given in ref. 7 (equations are in terms of frequency, 
v,rather than wavelength, X) To solve for i^ and the following approximations 

were made. (1) The temperature, Te, is uniform so that a^, and film index of 
refraction, n^, are also constant. (2) The film is a pure dielectric; therefore n f is a 
real number. (3) Boundaries at x = 0 and x = d behave diffusively; therefore i^(0), 

i x (d), reflectivity, p^., and emittance, e^. at the substrate-film interface and 

reflectivity at the upper boundary, are independent of 0. (4) Scattering is 

isotropic. (5) Interference effects are neglected. 

In solving the one dimensional radiative transfer equation (eq. (4)) it is convenient 

to split the intensity into two parts; i^ + for radiation in the positive x direction (0 < 0 

< Jt/2) and i^ - for radiation in the negative x direction (7t/2 < 0 < 7t). Using the 
definitions of i^ + and i^- and assumption (3) the boundary conditions at x = 0 and 
x = d are the following. 

x 

q£(0) = 2 ti ( 2 i£(0) cos0sin0d0 = 7ri£(0) = e^e^+p^q^O) 

Jo (7) 

qji(d) = ^x(d) = 7tpXoix(d) (8) 

Appearing in eq. (6) is the black body emissive power, 

27thCo 

SjlhC^T) = = — c 

X 5 [exp(hc 0 /XkT) - 1] /q% 


where h is Plank's constant, k is Boltzmann’s constant, and c D is the vacuum 

speed of light. 

Defining the optical depth. 



cc^dx' = a^x 


( 10 ) 


equation (4) can be integrated to obtain i^K^) and i^'(0), where 
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( 11 ) 


and those results together with the boundary conditions can be used to calculate the 
radiation fluxes(7) at the boundaries. 


qx(Kxd) = 2kJ* ix(Kxd,6)cosesined0 

= 2[£ Xs e Xb + P^€(0)]E 3 (K w ) + <X> + 
q^(0) = ~2nj % ix(O)cos0sin0d0 

=2pxoqx(Kxd)E 3 (Kxd) + «L 
Solving equations (12) and (13) results in the following 

+ rr n 2[ e Xs e w, + Pxs ^JE^Km) + 
qx(Kid) = ; ~ _2/t^ % 

1 - 4 PXsPXo E 3 ( K Xd) 


(12) 


(13) 


(14) 


qx(0) = 


2pXot^ + 2E 3 (K^)£ Xs e Xb ]E3(K Xd ) +3>_ 

1 “ 4 PXsPXx> E 3( K Xd) 


(15) 


where. 


f K xd * * * 

<b + = 2*1 S x (Kx)E2(K w - Ky)dK x 
<D_ = 2^ KW S x (Kl)E 2 (Kj)dKx 


(16) 


(17) 


and the exponential integral, E n (u), is 
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( 18 ) 


E n (u) = f v^expC- — )dv 

Jo V 


Once the source function, S^, is known, the boundary fluxes q^Kj^ and q^_'(0) 

can be calculated and thus the hemispherical spectral emittance as well. The 
emittance is defined as 


. _Qx(Kxd) 


(19) 


where Q^(K^j) is the spectral emissive power at = K^. To include the 

refractive properties of the boundary at x = d the maximum angle, 0 M , that allows 
radiation to escape the film is given by Snell's law. For x > d we assume n 0 = 1, 
therefore Snell's law yields the following result. 

= cos20 m = 1 " ( n x) 2 (70) 


For 0 > 0 M radiation is totally reflected. Including angles greater than 0 M would 
allow the possibility of > 1 (6). Therefore, using the boundary condition given 

by eq. (8) the emissive power, Q^, is the following 

/*®M 

Qx( K Xd) = 2 tc [ ix(K Xd )cos9sinede - p^ ii + (K^osOsinOde] 

Jo Jo (21) 

= (1 - PxJ(ql(Kxd) - + p fa qi(0)]E 3 (-H) . <^( 

r^M 

where 


/•K- 


^ - 271 V-m 


>Xd 


* Kj^-K* , 

Sx(K^)E 2 (— )dK x 

fm 


( 22 ) 


In reference 7 the source function is obtained for isotropic scattering and diffuse 
boundaries (assumptions (3) and (4)). The solution is obtained in the following 
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manner. First, the exponential integral functions in the source function equation are 
approximated by simple exponential functions. This allows the integral source 
function equation to be reduced to a linear second order differential equation which 
can then be solved. With no scattering the source function has the following simple 
form. 


n^ 2 Te) 0 


(23) 


Since we have assumed Tg is uniform, has no dependence and the 
integrations appearing in d>_ and can be carried out. Also, since the 
substrate is opaque 


Pxs~ 1 " E Xs 


(24) 


Therefore using equations (14), (15) and (24) in equation (21) yields the following 
result for the spectral emittance 


£ X = ^^{2h_[ £xs + n£(l - e^) (1-2E 3 (K w )] + ni[h + (l-2E 3 (K w )) - ^(l^E^))] 

(25) 


where. 


0 

h+ = 1 ~ 4 PXsPXol A M E 3( K Xd) E 3('^) 


h_ 


= E 3 (K Xd )-pJ I E 3 ( 


Kxd 

P-M 


) 


Dx=l-4pXsPxo E 3(Kxd) 


(26) 

(27) 

(28) 


At the emitting surface (x = d) we assume that the reflectance, p^, can be 

approximated by the normal reflectance for a specular reflective dielectric interface. 
Therefore, assuming n 0 = 1.0 for x > d the reflectance is the following(6) 


1 ^2 


(29) 
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Equation (25) shows the exponential dependence P^x) = 1/2 e(‘ 3/2x )] of on 
the optical depth, K^. Also, the importance of substrate emittance is apparent. For 

X regions outside the emission band where is small [E 3 (0) = 1/2] the emittance 

will be dominated by the substrate contribution. Thus it is important to have low 
substrate emittance to minimize the out of band emission and thus maximize the 

emitter efficiency. Low also implies p^. approaches 1 which means that a 

significant contribution (1 - e^. term in eq. (25)) to the emission band results from 

reflection at the substrate. Therefore, low substrate emittance minimizes out of 
band radiation and increases within band radiation. 



Optical Depth, K x = cx^d 

FIGURE 2. Thin film spectral emittance as a function of optical depth 
for scattering albedo=0, .5 and substrate emittance=.01,.2. Also, film 
index of refraction=1. 9, reflectance at film-substrate interfaced - 

In figure 2 the spectral emittance is shown as a function of the optical depth for 

two substrate emittances (e^ = .01, .2) and n^ = 1.9, which is representative of 

rare earth emitters. Results for both no scattering (cr^ = 0) and equal scattering and 

absorption coefficients (a^ = a^, = .5) are shown. The £>^=.5 results were 

obtained from reference 7. Several observations can be made from the results in 
fig. 2. First, scattering obviously reduces the spectral emittance however, emitter 
efficiency will not be greatly affected since the ratio of the emittance within the 
emission band to the emittance outside the emission band determines the emitter 
efficiency. This will be discussed later.. The second thing to notice is how the 
substrate emittance affects the spectral emittance. Although it is difficult to see in 
the figure, for K^, 0 the spectral emittance increases with increasing substrate 

emittance.However,as increases the opposite result occurs; the spectral 
emittance decreases as substrate emittance increases. This result occurs because of 
the reflectance term, (1 - e^), in the spectral emittance (eq. (25)). For 0 the 
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reflectance term is negligible. But as increases the reflectance term becomes 

important so that increasing reflectance (decreasing substrate emittance) results in 
increased spectral emittance. Thus, as already mentioned, low substrate emittance 
is important for two reasons. For emission outside the emission band where is 

small, low substrate emittance is required to minimize the emission. Within the 
emission band where K w is larger, low substrate emittance results in high 

reflectance at the substrate and thus increases the emission. 


EMITTER EFFICIENCY 


As stated earlier, the emission properties of the rare earth emitters can be 
approximated by a three band model and the emitter efficiency is defined as 


Qb 

116 Qi + Qb+Qo (30) 

and the power in each of the three bands is defined by 


f * 1 

Qb = ^Xb^ 1 


(31a) 

Qi= f e x e XbdX 
A 

(31b) 

Qu = J o e X e XbdX 

(31c) 


For the below and above band regions we assume a constant extinction coefficient, 
as well as, constant substrate emittance and film index of refraction so that the 
emittance is a constant in those bands. 


£1 = E x (K 1 ,£ ls ,n 1 ) 


Xj < X ^ °° , Kj — ctjd 


(32a) 


£ u — £ X (Ku’^us’^u) 


0 < X < X u , K u = a u d P2b) 


Here £i s and are the substrate emittances in the above and below band regions 
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nj and n u are the film indexes of refraction in those bands. For the emission band 
the extinction coefficient varies considerably for the rare-earth garnet emitters(5). 
Therefore, will be a function of wavelength within the emission band. Since £„ 

and £j are independent of X, the integrations in equations (31b) and (31c) can be 
performed so that the following is obtained for % 


t| e = U 


^Gfu!) + e u H(u„) , 
+ FO^X^Te) ] 


(33) 


where. 


G(u >>=fj 



dx 


hc ° 

Ul “X,kT E 


(34) 


H(u u ) = 

Ju 0 



dx 


hc D 

u ~^t e 


(35) 


r x i 


fo.,a u .te)=>.; 


e l(K>.> e Xs’ n Xi>) 


dX 

X 5 [exp(hc 0 /XkT E - 1] 



(36) 


The emitter efficiency for the thin film rare earth- YAG emitters was calculated using 
eq. (33). For no scattering the emittance is given by equation (25). With scattering 
the emittance was obtained from ref. 7. The extinction coefficients for Er-YAG and 
Ho- YAG were obtained from ref. 5. From the experimental extinction coefficient 

data the emission band is defined by X} = 1.66 pm (cq = .181 cm -1 ), X u = 1.42 pm 
(a u =.397 cm -1 ) for Er-YAG and by Xj = 2.15 pm (cq = .056 cm* 1 ) X u = 1.82 pm 
(a u =.61 1cm) for Ho- YAG. 
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FIGURE 3. Dependence of Er-YAG and Ho-YAG emitter efficiency, n E ,on scatterig 
albedo, G.and substrate emittance, £ t ,forT E =1500K. (index of refraction = 1.9) 

Figure 3 shows the efficiency for Er-YAG and Ho-YAG as a function of film 
thickness, d, at T E = 1500K for no scattering and for equal absorption and 

scattering coefficients (H = .5). The substrate emittance and scattering albedo were 
assumed to be constants (e s = £j s = £bs = = ^u)* Also, the 

refractive index was assumed to be constant (n^ = n^ = n^ = 1.9). The value n^ = 
1.9 is representative of rare earth oxides(9). Several observations can be made 
from the results shown in figure 3. First, there is an optimum film thickness for 

maximum efficiency. Second, scattering causes only a small reduction in r) E . 
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However, a thicker film is required to obtain maximum % when scattering exists. 
As will be discussed later, scattering has a more significant effect on the emitted 
power. Finally, to have large tj e the substrate must have low emittance. 

Now consider the effect of emitter temperature, T E , on the emitter efficiency. 
In figure 4 q E for Er-YAG and Ho-YAG is shown for T E = 1200K, 1500K, 
1750K and 2000K. The values for dj, a u , and n f are the same as in figure 3. 

For the scattering albedo, = fiy = 1/3 (o^ = 1/2 a^) was chosen as 

representative of the rare earth- YAG emitter. Figure 4 shows the strong temperature 
dependence of q E . Two observations can be made. First of all the efficiency 

increases rapidly with temperature until a maximum value is reached and then 
decreases. In the case of Ho-YAG maximum efficiency occurs for T E = 1600K. 

For Er-YAG maximum efficiency occurs for T E > 2000K. In references 7 and 10, 
where the emission band is characterized by a single constant extinction coefficient, 
maximum efficiency was found to occur for Eb/kT E = 4, where E b = hc 0 /X b . In the 

case of Er (X b = 1.54 nm) this corresponds to T E = 2400K and for Ho (X b = 1.99 
nm) Tg ** 1800K. Second, the thickness for maximum efficiency decreases as the 

temperature increases. For Er-YAG the effect is small. However, for Ho-YAG the 
change is significant At T E = 1200K the optimum thickness is approximately 0. 1 

cm but at T E = 2000K the optimum thickness is approximately 0.05 cm. The 
difference in the extinction coefficient within the emission band accounts for the 
different behavior between Er-YAG and Ho-YAG. The magnitude of a x for Er- 

YAG is approximately a factor of 2 larger than <x^ for Ho-YAG 5 within the 
emission bands. Also, there is a difference in the shape of the vs X curves for 

Er-YAG and Ho-YAG(5). Both materials have two peaks in the vs X curves. 

However, for Ho-YAG the peaks are separated by a greater wavelength difference. 
As a result the optimum thicknesses for Er-YAG are smaller than optimium 
thicknesses for Ho-YAG. 
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0 0.04 0.08 0.12 0.16 0.2 
Film Thickness, d cm 

4a) Er(40%)-YAG (X= 1.4^im, a u = .397 cm' 1 , 1.66nm, a= .181 cm’’) 



Film Thickness, d cm 

4b) Ho(25%) - YAG {\= 1.82, a = .611 cm' 1 , X= 2.15jim, a= .056 cm' 1 ) 
FIGURE 4. Emitter Efficiency of Er-YAG and Ho-YAG as a function of 
film thickness for various temperatures.Scattering albedo(Q=1/3), substrate 
emittance(e j =.1), and index of refraction(n=1 .9), the same in all bands. 
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5a) Er(40%) - YAG (X = 1.42*im, a = .397 cm* 1 , 1.66^m, a = .181 cm' 1 ) 



FIGURE 5. Power density as a function of film thickness for various emitter 
temperatures. Scattering albedo(n=1/3),substrate emittance(e j =.1), 

and index of refraction(n=1 .9) the same in all three bands. 


EMITTER POWER DENSITY 


Besides efficiency, the other important emitter performance parameter is the 
emitted power in the emission band, Q b . The efficiency is strongly dependent on 

the emittance ratios, and 8^/^, where e b is the emittance within the emission 

band. Whereas, efficiency is not strongly dependent on the magnitude of the 
emittance in the three bands (see eq. (33) and ref. 10). However, Q b is directly 
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proportional to the magnitude of the emittance within the emission band. 
Therefore, if the magnitude of the emission band emittance, £ 5 , is low, the emitted 

power, Q b , will be low. However, the efficiency can still be high if and £i/£b 

are small. The efficiency and power density are not strongly coupled. It is possible 
to have both high efficiency with low emitted power and low efficiency with high 
emitted power. 

In figure 5 the emitted power density, Q b , is shown as a function of film 
thickness for Er-YAG and Ho-YAG for the same conditions as figure 4. As can be 
seen, Q b increases monotonically with d. This results because £f, increases 

monotonically with d (see fig. 2). Thus by increasing thickness, emitted power 
increases, but efficiency will begin to decrease when the thickness exceeds the 

thickness for maximum efficiency. Also, since increasing scattering reduces (see 
fig. 2), Q b will also be reduced by increased scattering(7). Note that Q b for Ho- 
YAG is greater than Q b for Er-YAG for the temperatures (T E < 2000K) shown in 
figure 5. If the temperature were increased beyond 2000K, Q b for Er- YAG would 
eventually exceed Q b for Ho-YAG. 


CONCLUSION 

A radiative transfer analysis that includes isotropic scattering and assumes 
constant temperature throughout the film was used to calculate the spectral 

emittance, e^, of a thin film. Using based on measured extinction coefficients 
and assuming a three band model for the rare earth- YAG thin film emitter, the 
emitter efficiency, T[g, and emitted power density, Q b , were calculated. Efficiencies 

of approximately 50% are possible with power densities of several watts/cm 1 2 3 4 in the 
temperature range 1200 < T E < 2000K for Ho-YAG and Er- YAG. In order to 
attain high efficiency, low substrate emittance is required. Also, there is an 
optimum thickness in the range .04 < d ^ 0.1 cm to attain maximum efficiency. 
However, power density is a monotonically increasing function of film thickness, 
d. Finally, scattering has negligible effect on efficiency but a significant effect on 
power density. 
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